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I, Daniel K. Sodickson, do hereby declare that: 

1 . I am the applicant of the above-identified patent application and its priority 
application, and inventor of the subject matter described and claimed therein, 

2. Prior to tiie date of December 23, 1 999, while working iji the United States, I 
conceived the idea of a method and apparatus for implementing parallel MR image 
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3 . As recorded in Exhibit A, I had by a date prior to December 23, 1 999 already 
reduced to practice the invention set forth in the clainas of the above-identified application, 
which are all drawn to aspects of my invention set forth in the Exhibit A disclosure. 
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made are punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the 
United States Code, and that such willfixl false statements may jeopardize the validity of the 
applictaiion, any patent issuing thereon, or any patent to which this verified statement is 
directed. 
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1 . TITLE OF INVENTION : (Brief, but comprehensive, technically accurate and descriptive) 
A Generalized Basis Approach for Improved Parallel Magnetic Resonance Imaging 



2. IDENTIFICATION OF CONTRIBUTORfS) : Please provide the following information for 
each inventor: 



a. Full name 

b. Official title or position 



c. 
d. 

e, 
f. 

g- 



Department and Division 
Business phone number 
Fax number 
Email address 
Home address 



h. Citizenship 



Daniel K. Sodickson, MD, PhD 

Instructor in Medicine, Harvard Medical School 

Research Scientist, Beth Israel Deaconess Medical 

Center 

Department of Medicine, Cardiovascular Division 

(617) 667-8784 

(617) 975-5480 (FAX) 

dsodicks@caregroup.harvard.edu 

1 105 Massachusetts Avenue, Apt. 12-D 

Cambridge, MA 02138 

USA 



3. FINANCIAL SUPPORT/ CONTRACT IDENTIFICATION , This information is needed to 
determine whether this invention is subject to any conmiitments or restrictions arising from the terms of 
sponsorship. 

A. Identify the specific grant or contract number(s) (not the account number) and.the external sponsors 
(government agencies, industrial sponsors, private agencies or others) which provided support used to 
defray costs related to the research from which the invention resulted. 

• NTH Grant #R29-HL60802 

• Whitaker Foundation Biomedical Engineering Research Grant 

B. Was a biological, chemical or physical material or substance obtained from others to create this 
invention? 

YES NO X 

If yes, did a Material Transfer Agreement or other document accompany the transfer? 
YES NO 
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4. Please attach a CONCISE DESCRIPTION OF THE TECHNOLOGY , Your disclosure 
should enable someone having knowledge of the field to understand the invention and should cover ALL 
of the following points: 

a. General purpose and a technical description of what the invention is 

b. Describe the advantages and improvements over the existing technology 

c. Identify commercial applications 

d. Describe what is presentiy available in the field identifying existing methods, devices or materials (and 
their shortcomings). 

e. List published material such as patents, commercial literature and scientific articles relating to the 
invention. Search the internet for patents using: 
vyww.USPTO.gov (abstracts only) 
www.IBM.com (abstracts and claims) 

pctgazette.wipo.int. (foreign applications which are published 6 months after being 
filed) 

Introduction 

Imaging speed is of critical importance in many clinical applications of magnetic 
resonance imaging. Imaging of the cardiovascular system, for example, must contend with 
the challenges of cardiac motion, respiratory motion, and blood flow. Improvements in MR 
image acquisition speed over the past decade have allowed significant advances in the 
visualization and characterization of moving structures. Nevertheless, physical and 
physiologic constraints on the switching rate of magnetic field gradients and radiofrequency 
(RF) pulses place certain basic limits on the speed of traditional MR image acquisitions. 
Many modem MR scanners already operate close to these limits. 

Recently, a new approach to MR image acquisition has been developed, which can 
extend or circumvent the speed limits traditionally associated with gradient and RF hardware. 
This approach, known as "Parallel MRI," uses spatial information from arrays of RF coils toi^ 
perform some portion of the spatial encoding normally accomplished using gradients and RF 
pulses. The use of multiple effective detectors has been shown to multiply imaging speed, 
without increasing gradient switching rate or RF power deposition. 

This new paradigm of parallel MRI had its practical beginnings with the introduction 
of the SMASH technique (1-3), which allowed the first accelerated in vivo parallel images to 
be obtained. Subsequently, several related parallel imaging techniques have been proposed 
and applied to in vivo imaging. Parallel imaging techniques have tended to fall into one of 
two general categories, as exemplified by the SMASH (1-3) and the SENSE (4, 5) 
techniques, respectively. SMASH operates primarily in /:-space, while SENSE operates 
primarily in the image domain. Each approach has its particular advantages and limitations. 

The work described here represents a generalized formalism in which various 
techniques such as SMASH and SENSE appear as special cases. The generalized approach 
allows direct comparison of the two classes of techniques, and provides insight into origins 
of some of their distinct practical features. Appropriate modification of parameters in the 
generalized approach also allows the generation of hybrid techniques combining some of the 
advantages of the fc-space and the image domain approaches. Furthermore, numerical 



conditioning may be incorporated conveniently into the generalized formalism. This 
important feature of numerical conditioning enables effective compensation for noise or for 
systematic errors in coil sensitivity calibration, resulting in improved image quality Aind 
signal-to-noise ratio under a wide range of conditions. Incorporation of numerical 
conditioning directly into the parallel image reconstruction using the generalized approach 
also removes a cumbersome and potentially error-prone sensitivity calibration step involving 
division of two distinct reference images in the SENSE technique. 

In short, the generalized approach described here not only provides a unifying 
perspective connecting various parallel imaging techniques, but, in combination with 
appropriate numerical conditioning approaches, it also allows substantive improvemenis to 
techniques such as SMASH and SENSE, and enables hybrid approaches with further 
enhanced performance. 

Theory ^ 

We begin with the recognition that measured MR signal data is comprised of a 
generalized series of projections of the underlying distribution of MR-active spins in the 
imaged volume. A listing of the various projection functions then leads to practical strategies 
for parallel image reconstruction. 

"^The MR signal detected in any given RF coil is the result of a spatial integration of 
spin density (ignoring relaxation) against the sensitivity C/ of that coil and against the 
sinusoidal spatial modulations generated by encoding gradients: 

Si(K,ky) = lldxdypix, y) •Q(x,y)expHV)expHfc^y) (1) 

Here, / = 1,2,- is the index of any component coil in an L-coil array, and kx = 0,l,...,Nx-l 
and ky^OJ,.. .Ny-1 are thejk-space indices representing different frequency- or phase- 
encoding gradient steps. In other words, the signal comprises integrations or projections of ( 
the spin density against L*//x*iyy distinct functions ^' 

B,{x, y.k,,k^) = y)exp(-iV)expHfe^y) (2>:: 

The integral in Eq. (2) can be approximated with a discrete sum as follows, using a single k- 
space index k = (kx,ky) and a corresponding single pixel index; = (x,y): 

5,(fc„fc,) ^S^^J^Pj' Cj, cxp{-27Uijk) = j^pj' Bj^ (3) 

j i 

For any given coil /, this is just a discrete Fourier transform QWl), which is easily inverted 
with an inverse DFT: 

Notice, however, that if the signals from all the coils in the array are grouped together into a 
single index p = (k,l), a matrix equation is obtained for the spin density without modulation 
by any coil sensitivity: 

P,'-LS,B- (5) 
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In other words, if B can be calculated, the spin density can be determined. 

If a complete set of A^y phase encoding gradients are applied, the B matrix has 
dimension {Nx^Ny'^L) x {Nx'^Ny), This is clearly overdetermined. Phase encoding gradient 
steps may be omitted up to a maximum factor of L, and B will remain invertible in principle. 
In other words, when spatial information from an array of coils is available, the spin density 
may be determined from a reduced set of phase encoding gradients. This is the basis of 
spatial encoding in parallel imaging. 

. The structure of the B matrix is indicated graphically in Figure 1. In general, B will 
have dimension {Nx^Ny'^UM) x {Nx^Ny\ where M is the fraction of omitted gradient steps. 
The example in the figure shows a segment of the B matrix at a given x position, with 3 
component coils and a set of 8 phase encode gradient steps. Let us say that we wish to 
generate a reconstructed image with the equivalent of 16 phase encode gradient steps (i.e. we 
specify that M = 2). 

The rows of B in this example represent 3*8 = 24 vectors in a 16-dimensional space. 
Inversion operates by converting these vectors to a basis of 16 vectors spanning the space. 
Particular choices of the target basis result in inversion strategies which may be shown to be 
equivalent to SMASH or to SENSE. Altematively, hybrid approaches may be designed by 
choosing intermediate target bases, or by applying various approximations to the inversion. 




Figure 1: Graphical depiction of a B matrix, with different basis transformations leading 
to fc-space or image domain reconstructions, respectively. 



Image-domain reconstructions: the pixel basis 

Direct inversion of the B matrix yields a matrix which when multiplied by B produces 
the identity matrix. From the perspective of basis transformations, this represents 
transformation to a basis made up of distinct delta functions at each pixel of the image. For 
particular choices of inversion strategy, this may be shown to reduce to a SENSE image 
reconstruction, or to a reconstruction using the related subencoding (6) technique. 

k-space reconstructions: the Fourier basis 

Altematively, the Fourier transformation may be factored out of the B matrix, and the 
tranformed B matrix may be inverted independently from the FT matrix. This procedure 
corresponds to transformation to a Fourier basis made up orthogonal sinusoidal modulations. 
One further simplication yields a procedure equivalent to the SMASH image reconstruction. 



The transformed B matrix may be shown to be nearly block diagonal in many cases, as 
shown in Figure 2. SMASH image reconstruction corresponds to inverting only a subset of 
the transformed matrix, and assuming that the results apply equally well to all sub-blocks. ^ 
This assumption, which is based upon emulating the physics of phase encoding gradients, 
serves both to reduce image reconstruction time and to alleviate potential numerical 
instabilities which may result from inversion of the larger matrix in the presence of noise or 
systematic sensitivity calibration errors. 

Hybrid reconstructions 

One type of hybrid reconstruction results simply from expanding the inverted sub- 
blocks described above. Such an approach is shown schematically in Figure 3. These hybrid 
reconstructions share some of the numerical stability and efficiency of the SMASH 
technique, while also allowing some of the pixel-by-pixel control inherent in the image 
domain techniques. 

Another form of hybrid reconstruction results from using a target basis intermediate 
between the pixel basis and the Fourier basis described above. For example, a partially 
localized basis such as a wavelet basis may be used. 




Figure 2: Simulated modulation 
function matrix B for six-fold 
acceleration (M = 6) with a six-element 
coil array (L = 6). The Fourier 
transform has been factored out of B, 
resulting in a nearly block diagonal 
matrix. 
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I Figure 3: Schematic illustration of hybrid reconstruction 
I approach using inverted sub-blocks of varying size. Solid 
I boxes: minimum sub-block size corresponding to SM^^SK 
I reconstruction. Dotted box: maximum sub-block size ' 
I corresponding to inversion of the full B matrix in image 
domain reconstructions. Dashed boxes: expanded sub- 
i block size corresponding to hybrid reconstruction. 



Numerical conditioning 

In addition to the numerical conditioning provided by sub-block diagonalization, 
further numerical conditioning may be added directly to the B matrix inversion. 
In the overdetermined case, when B is rectangular with M < L, a generalized inverse 
procedure must be used. One such procedure involves performing a singular value 
decomposition (SVD) of ttie B matrix 

B = USV (6) 
where U and V are unitary matrices and S is a diagonal matrix containing eigenvalues of B in 
order of decreasing size. The inverse may then be performed as 

B-*=V'-S-'-U* ' (7) 

The smallest eigenvalues in S represent potential numerical instabilities in the ^ 
reconstmction. Since small eigenvalues will bfe inverted to large values in the inverse S* , 
either noise or systematic errors in sensitivity calibration affecting the eigenvectors 
associated with these small eigenvalues will be amplified in the reconstruction. Such 
amplification may be prevented by establishing a minimum eigenvalue threshold and 
performing one of the following procedures: 

1) Eliminate all eigenvalues below the eigenvalue threshold from die inversion. 

2) Set all sub-threshold eigenvalues equal to the threshold value. 

3) Add the threshold value to all eigenvalues. 

The result of this process is to modify or eliminate from the inversion the components most 
responsible for noise and error propagation, resulting in improved signal to noise ratio and 
reduced artifact in the reconstructed images. 

The price of these improvements is the omission of certain features of the 
reconstructed image in regions most susceptible to numerical instability. Unlike the 
sensitivity threshholding procedure described in (4, 5), which eliminates all regions falling 
below an established intensity and/or density threshold, the SVD approach results in an 
automatic threshholding strategy which only excludes or approximates regions of the image 
which are destined to yield large errors in the reconstmction. As a result, it may be used in 
conjunction with simple in vivo sensitivity references, without requiring additional 
processing or acquisition of an accompanying body coil image. 



Eigenvalue conditioning may be applied to B matrix inversions in the image-domain 
pixel basis, the k-space Fourier basis, or any mixed basis. Furthermore, if an alternative 
generalized inverse procedure is used, such as the least-squares pseudoinverse 

Bi«,=(B^Bf (8) 

eigenvalue conditioning may still be applied through an SVD of the inverse (b^b] . 
Results 

Figure 4 shows several examples of generalized basis parallel image reconstructions 
in a phantom imaging experiment. 




Figure 4: a) Reference image. b)-e) Six-fold accelerated phantom images in a double- 
oblique image plane, b) SENSE reconstruction, processed sensitivity reference (surface coil% - 
and body coil), c) Subencoding reconstruction, raw surface coil sensitivity reference, d) 
Generalized basis reconstruction with SVD eigenvalue threshold, raw surface coil reference, 
e) Generalized basis extended sub-block hybrid reconstruction, raw surface coil reference. 

Practical issues 

The dimensions of the B matrix may become quite large when images of high spatial 
resolution and large matrix size are acquired, and the inversions may consequently become 
both memory-intensive and time-consuming. Since the x and y integrals in Eq. (1) above are 
interchangeable, however, a separate smaller inversion may be performed for each position in 
the non-coil-encoded direction, and parallel image reconstruction may proceed line by line. 
This can represent a significant savings in image reconstruction time. 
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Areas of potential claims 

Method and apparatus to implement generalized basis transformations for parallel 
MR image reconstruction, which may be made to converge to SMASH or 
SENSE/subencoding, or to hybrid reconstructions which balance the advantages of each to 
yield higher quality images than would be possible with either alone. 

Numerical conditioning methods for the reconstructions described above, to improve 
noise and error propagation and to allow use of in vivo sensitivity references without 
thresholding or extensive post-processing. 

Tailored image reconstmction hardware and software, which may be combined with 
tailored RF coil arrays and interfaced with clinical and research MRI systems for improved 
imaging speed and efficiency. 

Additional information 

For additional information, please see the attached abstract (p, 10), which has been 
submitted for an upcoming conference of the International Society for Magnetic Resonance 
in Medicine in April of 2000. 



References 

1) Sodickson DK, Manning WJ, Simultaneous acquisition of spatial harmonics (SMASH): 
fast imaging with radiofrequency coil arrays. Magn Reson Med 1997;38(4):59 1-603. 

2) Sodickson DK, Griswold MA, Jakob PM. SMASH Imaging. Magn Reson Imaging Clin 
NAml999;7(2):l-18. 

3) US Patent #5,910,728, June 8 1999, Sodickson DK inventor, BIDMC assignee. 
Simultaneous Acquisition of Spatial Harmonics (SMASH): Ultra-Fast Imaging with 
Radiofrequency Coil Arrays. 

4) Praessmann KP, Weiger M, Scheidegger MB, Boesiger P. SENSE: Sensitivity encoding >^;> 
for fast MRL Magn Reson Med 1999;42(5):952-962. 

5) International patent application #PCT/IB99/00652, international publication 
#WO 99/54746, October 28, 1999. Pruessmann KP, Weiger M, Scheidegger MB, 
Boesiger P inventors. Philips Electronics assignee. 

6) Ra JB, Rim CY. Fast imaging using subencoding data sets from multiple detectors. Magn 
Reson Med 1993;30(1): 142-145. 



9 



A Generalize^ Jasis Approach to Spatial Encoding Coil Arrays: 
SMASH-SENSE Hybrids and improved Parallel MRI at High Accelerations 

D. K. Sodickson 

Department of Medicine, Cardiovascular Division, BeOi Israel Deaconess Medical Center and Harvard Medical School, Boston, Massachusetts. USA 



introduction 

The past few years have seen the emergence of a new paradigm of rapid 
MR imaging in which multi-component detector arrays are used to perform 
some fraction of image data acquisition in parallel. Parallel imaging tech- 
niques have tended to fall into one of two general categories, as exemplified 
by the SMASH (1) and the SENSE (2) techniques, respectively. SMASH 
operates primarily in Jt-space, while SENSE and the related "subcncoding" 
reconstruction technique (3) operate primarily in the image domain. Each 
approach has its particular advantages and limitations. 

Here a generalized formalism is presented directly demonstradng the 
precise relationship between SMASH and the image domain reconstructions. 
Not only does this formalism lend insight into the similarities and difTecences 
between the techniques, but it also allows the generation of hybrid recon- 
structions exploiting the advantages of both. In addition, it facilitates the 
introduction of numerical conditioning which may be used to improve the 
quality of parallel image reconstructions, especially at high acceleration fac- 
tors. In this work, elements of the generalized formalism are described, and 
its use is demonstrated in imaging experiments. 

Theory 

Each acquired MR signal data point may be viewed as a spatial integra- 
tion of the true image against the joint modulation functions imposed by 
applied gradients and component coil sensitivities: 

5;(k) = Jd V(x)C;(x)expH'k • x) = Jd'xp(x)B,(k,x) (1) 

Here / is the array element number, Cj is the sensitivity of coil /, p is the ap- 
parent spin density including relaxation effects, and k is the wavenumber of a 
periodic function resulting from evolution in frequency* and phase-encoding 
gradients. The general theory of sensitivity encoding outlined in (2) begins 
with a similar starting point. 

If the spatial integral in Eq. (1) is approximated with a discrete sum, we 
may write die signal Sn as 

or, in matrix form, 

S«Bp (3) 
The central step in a parallel image reconstruction now involves inverting the 
matrix B to recover the spin density p. 

In the discrete model of Eq. (3), the various modulation functions present 
in the rows of B may be recognized as vectors in a vector space which also 
contains the reconstructed image. In particular, for an acceleration factor of 
M, if (N/M) encoding gradient steps are used with an array of L component 
coils, and the desired reconstructed image is to be made up of pixels, B 
comprises a set of {N/M)*L vectors spanning an ^-dimensional vector space. 
In odier words, the acquired MR signal simply represents (N/M)*L different 
projections of the "true" image onto a mixed Fourier and sensitivity basis 
representing the target image. Various basis transformations may then be 
performed to facilitate and condition the inversion of B. 

If B is transformed to the Fourier basis, i.e. if the Fourier transform is 
factored out of the inversion, the relationship between Jk-space and image 
domain approaches becomes readily apparent. It may be shown that SMASH 
represents inversion of a small subset of the larger matrix. This inverted 
subset is dien replicated in a block diagonal manner to form die fiill inverse. 
SENSE, on the other hand, involves a particular approach to inversion of the 
full matrix with elements that have been smoothed, filtered and thiesholded 
using the sensitivity mapping procedure described in (2). 

In many cases, the transformed matrix is already nearly block diagonal, 
so that the SMASH reconstruction represents a compact and computationally 
efficient solution to die general problem. Figure 1 shows a contour plot of 
the eleinents of a simulated B matrix for such a situation. In cases for which 
the matrix does not have this property, the block diagonalization inherent in 
SMASH reconstruction provides numerical conditionmg at the expense of 
global accuracy. Rather than smoothing coil sensitivities, SMASH may be 
shown to entail a smoothing of the spatial variations among pixel-by-pixel 
weighting factors in the reconstruction, which serves to control the localized 
effects of noise and regional sensitivity calibration errors. 

Using the generalized basis approach described above, hybrid techniques 
may be formed by simply extending the size of the inverted sub-block. This 
allows a smoodi trade-off between the numerical conditioning mherent in the 
SMASH reconstruction and the pixel-by-pixel control of SENSE and suben- 
coding. 



Conditioning may also be inut)duced explicitly into the reconstruction 
using additional straightforward basis transformations. For example, a sin- 
gular value decomposition (SVD) of the B matrix yields the eigenfunctions 
of B in order of decreasing eigenvalue. Eigenvalue thresholds may then be 
applied to exclude from the inversion those componerts most responsible for 
noise and error propagation. Unlike the sensitivity thrcshholding procedure 
described in (2), which eliminates all regions falling bf'low an established 
intensity and/or density threshold, the SVD approach results in an automatic 
thrcshholding strategy v;hich only excludes or approximates regions of the 
image which are destined to yield large errors in the reconstruction. As a 
result, it may be used in conjunction with simple in vivo sensitivity refer- 
ences, without requiring additional processing or acquisition of an accompa- 
nying body coil image. 

Methods and Results 

Geneialized basis reconstructions were carried out and compared with 
pure SMASH or pure SENSE or subencoding reconstructions for phantoms 
and for in vivo cardiac images in a healthy adult volunteer. By appropriate 
tuning of parameters in the generalized approach, images were made pixel- 
by-pixel identical with eitiier SMASH or SENSE reconstructions. Alterna- 
tively, additional conditioning was added through hybrid reconstructions 
and/or eigenvalue thresholds. Figure 2 shows several double-oblique phan- 
tom image reconstructions at a fector of 6 acceleration - the theoretical 
maximum for die six-element coil array used in diese experiments. 

Discussion 

Numerical conditioning becomes particularly important at high accelera- 
tion factors, when the matrix inversions involved in any parallel reconstruc- 
tion strategy can most easily become singular. It is also important for large 
angulations of the image plane widi respect to die coil array, or in general for 
any situation in which die spatial distinctness of component coil sensitivities 
is degraded. The generalized basis approach ouUined here allowed improved 
reconstructions of double-oblique images at a high acceleration factor. It 
also allowed use of simple in vivo sensitivity references, preserving small 
regions of signal surrounded by noise, such as pulmonary blood vessels in 
thoracic images. Future work will involve optimization of the generalized 
basis approach for clinical applications of psuallel MRI at high accelerations. 
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Figure 1: Simulated modulation Junction matrix B far six-fold acceleration 
(M ~ 6) with a six-element coil array (Ls^6}. 




Figure 2: a) Reference image, b)-e) Six-fold accelerated phantom images 
in a double-oblique image plane, b) SENSE reconstruction, processed 
sensitivity reference (surface coil and body coil), c) Subencoding recon- 
struction, raw surface coil sensitivity reference, d) Generalized basis re- 
construction with SVD eigenvalue threshold, raw surface coil reference, 
e) Generalized basis extended sub-block hybrid reconstruction, raw surface 
coil reference. 
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5, PUBLIC DISCLOSURE / PUBLICATION PLANS : In the United States, a patent 
application must be filed no later than one year after public disclosure of the invention in 
detail. In other countries, filing must take place before public disclosure; however, where 
there has been a U.S. filing before publication, a one-year grace period is granted for 
foreign filing. Public disclosure includes abstracts and presentations at scientific meetings (including 
poster sessions), public seminars, publications, disclosure to others outside of the Medical Center who have 
not signed a confidentiality agreement, and use, sale or offer of sale of the invention. Identify dates and 
circumstances of any such disclosures. Also, indicate your future disclosure or publication plans, and 
NOTIFY the Office of Corporate Research if the invention becomes publicly disclosed or published in the 
future (whether by plan or inadvertently). 



The attached abstract has been submitted to the confidential review board for an 
international conference on Magnetic Resonance in Medicine to be held on April 1-7 2000. 
Abstracts for this conference are expected to be made public some time in March. 



6. Please attach an ABSTRACT OF THE INVENTION: This will be used in our campaign to 
promote the invention and should, therefore, be no longer than a typewritten page and not contain ^^ 
confidential information. :^ 

Please see the attached submitted conference abstract. 
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1. POTENTIAL LICENSEES . Provide as much detail as possible. 
a» List any commercial entities that may be interested in licensing this invention. 

List commercial entities, if any, that you specifically do not want contacted regarding this technology 
and please indicate why. 



1) Large Magnetic Resonance Imaging system manufacturers: Siemens, Philips, General 

Electric 

2) Small MR or radiofrequency coil companies 

3) Venture capitalists 



I (we) hereby agree to assign all right, title, and interest to this invention to Beth Israel Deaconess Medical 
Center (Medical Center) and agree to execute all documents as requested, assigning to Medical Center our 
rights in any patent application filed on this invention, and to cooperate with the Medical Center Office of 
Corporate Research in the protection of this invention. Medical Center will share royalty income derive^d 
from the invention with the inventor(s) according to its standard policies. 
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